Non-linear emission spectra of quantum dots strongly coupled to photonic mode 
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A theory of optical emission of quantum dot arrays in quantum microcavities is developed. The 
regime of the strong coupling between the quantum dots and photonic mode of the cavity is con- 
sidered. The quantum dots are modeled as two-level systems. In the low pumping (linear) regime 
the emission spectra are mainly determined by the superradiant mode where the effective dipoles 
of the dots oscillate in phase. In the non-linear regime the superradiant mode is destroyed and 
the emission spectra are sensitive to the parity of quantum dot number. Further increase of the 
pumping results in the line width narrowing being an evidence of the lasing regime. 

PACS numbers: 42.50.Ct, 42.50.Pq, 78.66.-m, 78.67.Hc 



I. INTRODUCTION 

The light-matter coupling on the nanoscale is in focus 
of the research during the last decade owing to the possi- 
bilities to realize quantum electrodynamical effects in the 
solid state and to the perspectives of the efficient emission 
control needed for device applications. Semiconductor mi- 
crocavities with embedded quantum dots are among the 
most promising objects in these respects. Two regimes of 
the light-matter interaction in such systems can be identi- 
fied. The first one, weak coupling regime, is characterized 
■by the relatively small interaction constant as compared 
with the decay rates of the photonic and excitonic states. 
In this case, the microcavity merely modifies the quantum 
'dot emission rate owing to Purccll effect. The second, 
strong coupling regime corresponds to the large coupling 
constant as compared with the decay rates which makes 
possible the coherent energy transfer between the quan- 
tum dot and a cavity resulting in formation of the mixed, 
half-light - half-matter modes called exciton-polaritons. 
Such a strong coupling regime was observed recently in 
semiconductor microcavities with quantum dots*^ 

Although, in most of the experimental and theoreti- 
cal studies the simplest possible situation of one quantum 
dot being close in energy to the cavity resonance was con- 
sidered so farj^ji it is feasible to achieve the strong cou- 
pling for multiple quantum dots£ Theoretical analysis^ 
demonstrated that, at the weak pumping where both cav- 
ity mode and quantum dots can be treated as classical 
harmonic oscillators the optical properties are dominated 
"by the superradiant mode which can be considered as 
a collective in-phase oscillation of quantum dot dipole 
moments^ In perfect systems this mode is the only one 
which interacts with a photon and their coupling constant 
is strongly enhanced as compared with the single dot case. 
The superradiant mode is relatively stable to the impact 
of disorder (i.e. spread of quantum dot resonant ener- 
gies) and can be even stabilized by the allowance for the 
interdot tunneling^ 

The increase of the pumping rate yields the increase 
of the electric field in a cavity and of the exciton dipole 
moment oscillation amplitudes. Thereby, the model of 
the classical oscillators becomes violated and the unhar- 
monicity effects have to be taken into account. Important 
questions in this respect are (i) how the emission spectra 
of the strongly-coupled system of quantum dots and cavity 



change as a function of pumping and (ii) whether the su- 
perradiant mode survives to some extent in the non-linear 
regime? The present paper is devoted to the theoretical 
study of these questions. 

We develop here a theory of non-linear emission of 
quantum dot ensembles embedded in microcavities in the 
strong coupling regime. We assume that the quantum dots 
are relatively small so that the electron and hole forming 
a zero-dimensional exciton are quantized independently 
and no more than one electron-hole pair can occupy the 
dot. The incoherent pumping of the excitonic states is 
considered. We also take into account the effects of the 
quantum dot resonant energies spread. 

The main results of our paper can be summarized as 
follows: 

1. The emission spectra of the single quantum dot de- 
tuned from the cavity resonance strongly depends 
on the pumping rate: it transforms from the dou- 
blet with dominant dot emission line to the single 
peak spectrum positioned near the cavity resonance 
with an increase of the pumping power. 

2. In the case where the spread of resonance energies of 
quantum dots is negligible, the non-linear emission 
spectra at moderate pumping rates exhibit either a 
single or a double peak structure depending whether 
the number of dots in a cavity is even or odd, respec- 
tively. 

3. At relatively high pumping rates the emission shows 
a single line whose width decreases with an in- 
crease of the pumping being a signature of the lasing 
regime. 

The paper is organized as follows: Section [Til describes 
the model assumptions and outlines the calculation pro- 
cedure. We demonstrate, that in strong coupling regime 
with well defined polariton states, the optical spectra can 
be found from the kinetic equation. An advantage of this 
approach over the full density-matrix calculations^ be- 
sides lower computational costs, is the possibility to de- 
rive transparent analytic answers for the optical spectra. 
In Sec. IIIII the developed method is applied to the case 
of a single quantum dot where the results known in the 
literature are reproduced in a simple fashion in the case 
of a dot being degenerate with the cavity mode. We also 
consider in this Section the situation of a dot detuned 



from the cavity mode and demonstrate the specific prop- 
erties of the emission spectra. Section lTVl is devoted to the 
case of several quantum dots in a cavity and the quantum 
dot number parity effects arc studied there in detail. To 
a certain extent, our results are valid in the case of a 
single quantum dot with multiple size-quantization levels 
in the vicinity of the photonic mode resonance, studied 
in Ref. |T(J. Concluding remarks are presented in Sec. |V] 
Derivation of analytical results for the microcavity with 
single dot is presented in Appendix [A] The effects of 
cxciton-cxciton interactions arc briefly discussed in Ap- 
pendix [B] 



II. MODEL 

Wc are interested in the microcavity with N embedded 
quantum dots. It is assumed that only one photonic mode 
of the cavity is of importance, i.e. it is close enough to 
the quantum dot transition energies. Each quantum dot 
is described by a two-level system where the ground state 
corresponds to the empty quantum dot and the excited 
state corresponds to the dot occupied by a single electron- 
hole pair (exciton) . Such an assumption is readily realized 
for the small quantum dots whose effective size is smaller 
as compared with the exciton Bohr radius ^ For the pur- 
poses of the present paper the spin degrees of freedom of 
photons and excitons are disregarded. 

Under above assumptions the Hamiltonian of the stud- 
ied system can be written in the following forroi^ 



A ! 



U = huje^c + Y^ Twx,ibl b i + h S^2( cJ(b i + cb 



(1) 



i = l 



Here loq is the resonance frequency of the cavity u>x,i 
(i — 1,...,N) are quantum dot resonance frequencies, 
c^, c are the creation and annihilation operators for the 
photon mode, respectively, and b], b, t are the analogous 
operators for the quantum dot modes. Cavity mode is 
bosonic while quantum dots are treated as two-level sys- 
tems. Therefore operators b\, bj obey Fermi commutation 
rules for i = j {b\b i + b^] = 1), while for different i =^ j 
quantum dot operators simply commuted The last term 
in Eq. (Q} describes the coupling between the excitons 
and photon, Kg is the coupling strength. It is assumed 
here that the coupling constant is identical for all dots, 
the developed theory can be easily extended to allow for 
different coupling strengths. 

The different nature of photon and exciton states makes 
it impossible to consider the dynamics of the coupled 
quantum dots/microcavity system in a model of classical 
harmonic oscillators. Indeed, the quantum dot cannot ac- 
commodate more than one exciton, although the number 
of photons in the cavity is not limited. At small pump- 
ing densities the excitation is transferred back and forth 
between the exciton and photon modes. An increase of 
the pumping power can result in the appearance of extra 
photons in the cavity. These extra photons cannot be ab- 
sorbed by the quantum dot but modify the light-matter 
coupling instead. 

In our model we neglect completely the interaction be- 
tween excitons. Clearly since wc model each dot as a two 



level system the number of excitons in each quantum dot 
is no more than one and there is no interaction within the 
same dot. The Coulomb interaction of excitons in different 
dots is neglected as well assuming that it is much smaller 
than the coupling constant Hg. The effects of interactions 
are briefly discussed in Appendix iBl 

The Hamiltonian (Q]) describes the eigenstates and the 
energy levels of the quantum dots in the microcavity in 
the absence of pumping and the non-radiative decay of 
the excitons and photons. Wc introduce the pumping rate 
to the ith exciton state, Wi, that is the number of exci- 
tons generated in a quantum dot per unit of time and 
the decay rates of the exciton and photon populations, 
Tx,i (i = l,...,iV), Tc, respectively. The cavity decay 
rate Tq defines the number of photons leaving the cav- 
ity per unit of time due to the non-zero transparency 
of the mirrors and the exciton population decay rates 
Tx,i characterize the non-radiative processes in quantum 
dots. In the present paper we focus on the strong cou- 
pling regime where the eigenstates of the Hamiltonian ([l) 
are well defined. Therefore, we assume that rx,;,rc <C g 
(i = 1 , . . . , N) . The requirement of the well defined states 
also imposes a certain restriction on the pumping rate 
W< g 2 /T c ,g 2 /T x ,i, see below. 

The Hamiltonian (|TJ) describes the energy transfer be- 
tween the quantum dots and the microcavity and con- 
serves the total number of particles in the system (excitons 
and photons). Therefore, its eigenstates can be labeled by 
the total number of particles, to, (and, of course, by other 
quantum numbers which take into account the degener- 
acy of non-interacting to particle states, e.g. to photons, 
excitons or to— 1 photons, 1 exciton, etc.) For instance, 
for a single quantum dot in the microcavity tuned exactly 
to the photon mode, ujx = ^>c the states are combinations 



(|m,0)±|m-l,l))/v / 2, 



(2) 



where |nc, nx) denotes the state with the definite number 
of photons, nc, and excitons, nx, with energies 



E rHi ± = mfiwc ± \fmhg. 



(3) 



These states with different m form Jayncs-Cummings 
ladder M^ The extension of Eqs. ©, © for the impor- 
tant cases of a detuned quantum dot and multiple quan- 
tum dots in a cavity are presented below in Sees. Mil and 
IVl respectively. It is seen from Eq. ([3]) that the energy 
spectrum of the interacting quantum dot and microcavity 
depends strongly on the number of particles. The light- 
matter coupling strength increases with to. 

The manifolds, i.e. the sets of states with fixed to, are 
intermixed by the dissipative processes. Namely, the emis- 
sion of photon results in the transition from the manifold 
with to particles to the manifold with to— 1 particle, while 
the generation of an exciton in a quantum dot leads to the 
increase of the number of particles by 1, i.e. to the tran- 
sition from the manifold to to the manifold to + 1. In the 
strong coupling regime one can describe the populations 
of the states in the framework of the distribution function 
fm,n, where the subscript \i enumerates different states 
in the manifold [e.g., u = ± in the case of single dot in 
a cavity, see Eqs. ©, ((3J)]. The distribution function is 



governed by the phenomenological kinetic equation which 
has the following form 



of the complete density matrix of the system Q m ,(i;m' ,(i> ■ 
It satisfies the following equation: 



d/, 



m,jU 



at 



rn' — m±l,/i' 



where D m ^ is the total decay rate of the state to, \x caused 
by the transitions to lower and upper manifolds resulting 
from the photon escape from the cavity, non-radiative ex- 
citon decay and pumping, 

m' — mil,/i' 

and W m ' ifi '^ m ^ are the corresponding transition prob- 
abilities. They can be separated into two parts, corre- 
sponding to the photon escape from the cavity, W c , and 
exciton generation in ith. dot, W^' 1 : 



W, 



m,fi^t-m r ,/i' 



ir 



m,fi— >rn- 






/i— >m' ,/i' ' 



(5) 
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n~, 



m,/j— >?n',/j' — ■*• COto'jTO— l|Cm— l,ji';m,;i| j 

- (r x ,^ m ', m -i + WiS m , tm+1 )\b^ im J 



rn,fj,—^m' ,fj,' 



l(<) 



and 6^/ ,. m arc the matrix elements of 



xiere Cm'^';^^ 

annihilation operators of photon and exciton in i-th dot 
respectively, taken between the states m, /x and to', /i'. It 
is worth to stress, that both pumping and decay events 
change the number of particles by one, thus coupling only 
neighboring manifolds: to, to ± 1. Note, that the man- 
ifold number to in Eq. (jj]) varies from to infinity, i.e. 
including the ground state of the system being the state 
with no particles. The distribution function obeys the 
normalization condition 



J2Ef> 



7n,fj, 



= 1. 



(G) 



m=0 /^ 



In the kinetic approach the emission spectra can be cal- 
culated by the Fermi golden rule 



I(u>) oc 



m,fi,fi' 



/, 






-l,n'&\Em,n — E m -i : 



Mi 

(7) 

and are determined by the rate of the photon escape from 
the cavity W^ M _> TO _ ljM ,. Here E m ^ is the energy of the 
corresponding state, and the (5-function ensures the en- 
ergy conservation. The common factor in Eq. ([7]) contain- 
ing the mirror transmission coefficients, etc., is ignored. 
To derive Eq. ([7]) we considered the process of the pho- 
ton escape from the cavity and represented the emission 
intensity as a weighted with the stationary distribution 
function / m .^ sum over all manifolds. Within the appli- 
cability of the kinetic equation the line broadening is to 
be small as compared with the transition energy, the con- 
dition surely satisfied in the strong coupling regime. It 
is seen that Eq. ([7]) gives infinitely narrow emission lines 
due to the presence of the <5-functions. In order to find 
the linewidths one has to go beyond the kinetic equation 
approach and take into account the non-diagonal elements 



-77 = fc^'^] + -^( 2c £ ct -JcQ- QC<c) + 

J2^(2b t gbl-btb t g-gbtb t ) + 

i 

^^■Wt&i-WU-QbA)- (8) 

i 

Here the first term describes the Hamiltonian-driven dy- 
namics of the coupled quantum dots and cavity system, 
second term describes the photon escape through the cav- 
ity mirrors, last two terms stand for the non-radiative 
exciton decay and non-resonant exciton pumping, respec- 
tively We use the general expression for the emission 
spectrum 16 



/(w)cxRc/ dte lut (c T (0)c(i)) 



(9) 



where c(t) and c'(t) are the time-dependent field oper- 
ators for the cavity mode, the angular brackets denote 
both quantum mechanical and statistical averaging and 
express the correlation function of the photonic mode via 
the density matrix as 



I(uj) CX ^ /l 



m,A» c m— 1,/j' ;m,fj, 



E 



^-?ni — l,/j,^;?ni ,/j-i • x - 



m,fi,fi' 






C™-i,;H- (io) 



Here 



~(m,/x,//) / ,\ _ p„ / A+Jut n (to,/*,/*') /.s 

Qmu^m-,-1^' Ke / dfe 0mi,pi;mi-lA.iW> 

where o„ ''*' „ , , (t) are the time-dependent solutions 
of the total density matrix equation ([8]) . The upper sub- 
scripts m, fi, fi' in Eq. (fTO]) label the manifold, populated 
at t — 0, so the initial conditions read 



(m, /*,//) 



-mi^iimi-l/i^ 



i(0) 



^77i,mi^/i,/ii 0p/,/^ 



(11) 



Thus, each term in the sum (|10[) with given indices to 
and fi is proportional to the stationary population f m ,fj,, 
similarly to Eq. (UJ. In derivation of Eq. (fTU)) we took 
into account that, in the strong coupling limit, the sta- 
tionary density matrix g st is diagonal in the basis of the 
manifolds, and can be found from the kinetic equation ((3]), 
o st , = S iS if 

Eq. (fTU|) is the particular case of the more general ex- 
pression presented in Ref. |9| which is valid for the arbi- 
trary coupling strength. In our case, the strong coupling 
condition results in the substantial reduction of the com- 
putational costs since much smaller number of the non- 
diagonal elements of the density matrix are to be cal- 
culated. If one neglects dissipative processes completely, 
then 



L> 



(ro,/*,/*') 

m-i ,/*i;mi — 1/*^ 



_ e -i(£»,„-B m -i, 



)'/ fi (5 8 S , 1 

u m,m± u fj,.fj,i v/j,' ifx^ , 



and Eq. dTDJ reduces to Eq. Q. 



III. SINGLE QUANTUM DOT 

In this section we analyze emission spectra of a micro- 
cavity with single quantum dot. Our calculation technique 
is tested in well-studied case of zero detuning between ex- 
citon and the photon mode^ ' 17 i 18 {loq = Wx) and then 
applied to a general case, wc ¥" w x- 
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FIG. 1: Photoluminescence spectra for N — 1 quantum dot 
in resonance with the cavity mode for different pumping rates 
W. Panels (a)-(e) correspond to W/Tc = 0.01, 1, 2, 30, respec- 
tively, and are calculated for o;x = ujc,^c/g = 0.1,rx/<? = 
0.01. Energies of the transitions, mainly determining the spec- 
tra, are indicated in panels (a), (b) by vertical lines. Dotted 
curve in panel (a) shows the Lorentzian with half-width equal 
to Fc/2, corresponding to empty cavity spectrum in the linear 
regime. Photon mode energy is used as the reference point, 
the spectra are normalized to their maximum values. 



A. Resonant quantum dot 

Photoluminescence spectra for zero detuning, calcu- 
lated for different rates of incoherent pumping rate W, 
are shown in Fig.[TJ Variation of the ratio between pump- 
ing and decay rates controls the population of the Jaynes- 
Cummings states, Eq. ^. Radiative transitions between 
these states determine optical spectra. Several regimes of 
light emission are resolved for different pumping rates: 

(a) Linear regime, W <C Tq, see Fig. QJi. The mean 
number of polaritons (to) ~ W/Tq is much less than unity, 
so that only the manifolds with to = and m = 1 are 
populated, /o 3> /i,±- Photon emission from two states 
with to = 1 results in the Rabi doublet in spectrum with 
peaks at 

Ei t ± — Eq = hiOQ ± hg . 

In these case photon and exciton modes can be treated 
as coupled harmonic oscillators since average number of 



photons and excitons in the system is very small and the 
statistics of these states becomes unimportant. 

(b) Intermediate regime, W ~ Tq- The mean number 
of polaritons is on order of unity and several lowest mani- 
folds are populated. The emission energies correspond to 
the difference of level energies in all relevant neighboring 
manifolds, hence, the spectrum has a complex multipeak 
stucture. In particular, for the case of Fig. QTb), where 
W = Tq, the peaks correspond to emission from man- 
ifolds with to = 1 ... 3, with the energies indicated by 
vertical lines. 

(c) High pumping regime, W ^S> Tq- The mean number 
of photons is much larger than one, (to) w W/(2Tc) 3> 1. 
Here, only the optical transitions between the states 

Eq. @ with the same symmetry, i.e. m,-\ > rn — 1, +, 

are important. The corresponding photon energies are 



E m ,± — £m-i,± ~ fiwc ± 



fig 



(12) 



resulting in a two-peak emission spectrum, see Fig. [TJ;. 
Simple analytical expression 



I {to) ex 



1 



(w - W - .9cff) 2 + 7 2 (w - loq + g cS ) 2 + 7 2 ' 

(13) 
well describes the spectrum, sec Appendix [S] for details 
of derivation. The doublet (fl~3|) is characterized by ef- 
fective Rabi splitting 2g e s = g/y/(m) and peak half- 
width 7 = {W + Tc)/4. Hereafter we neglect exciton 
non-radiative decay as compared with the photon leakage 
through the mirrors. We note, that although Eq. (|L3"|) 
has a simple form, it is determined by a lot of optical 
transitions Eq. (|12l) with different to. As shown in Ap- 
pendix [A] the interplay between these transitions due to 
the fermionic nature of the exciton leads to the broadening 
of the peaks with respect to the linear regime, cf. Figs.QJi 
andQJ;. Moreover, the distance between the peaks is much 
smaller than in linear case, because g e g <SC g. This effect 
allows semiclassical interpretation, since the mean number 
of cavity photons is much larger than unity. The decrease 
of <?eff with pumping can be considered as quenching of the 
optical transition due to the saturation of two-level dotJ^ 
In agreement with this concept, the average number of 
excitons is close to 1/2, just like for a two-level system, 
interacting with strong (classical) electromagnetic field^ 
(d) The further increase of pumping rate leads to the 
the transition to lasing regime, which takes place when W 
is equal to 



W* 



2g^T^ 



(14) 



The lasing regime is manifested as single-peak spectrum, 
with half-width rapidly decreasing when pumping rate 
grows. The approximate analytical expression for spec- 
trum reads 



-T(w) oc 



1 



(w - c^c) 2 + {g 2 r c /w*y 



(15) 



In agreement with this result Fig. [TJi demonstrates, that 
the emission spectrum becomes much narrower than in 
the linear case, cf. solid and dashed curves. The spectral 



width Tc(g/W) 2 is also smaller than the effective Rabi 
splitting 2g cS . 

Further increase of the pumping results in the self- 
quenching of a single-dot laser regime^^i and, corre- 
spondingly, violates the strong coupling condition. In- 
deed, if the pumping rate exceeds the splitting between 
the + and — states in the 'actual' manifold, W > y/(m)g 
or 



W > 



g 

rv 



the states + and — within the given manifold become 
ill-defined. It is worth noting that, in high finesse micro- 
cavities, with g 3> Tc, the characteristic pumping rates 
separating different spectral regimes are considerably dif- 
ferent, r c < W* < 5 2 /r c . Provided that W < g 2 /T c 
the strong coupling regime is maintained and our results 
are in the perfect agreement with those obtained in Ref. |9| 
by keeping all elements of the density matrix which pro- 
vides the justification of calculation technique presented 
in Sec. M 
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FIG. 2: Photoluminescence spectra for N = 1 quantum dot de- 
tuned from the cavity resonance for different pumping rates W. 
Panels (a)-(e) correspond to W/Tc = 0.02, 2, 10, 30, respec- 
tively, and are calculated for detuning cjx — Wc = 5g. Other 
parameters are the same as for Fig. [TJ Energies of the transi- 
tions, mainly determining the spectra, are indicated in panels 
(a)-(c) by vertical lines. Photon mode energy is used as the 
reference point, the spectra are normalized to their maximum 
values. 



B. Detuned quantum dot 

Now we turn to the case when the exciton energy is de- 
tuned from the photon mode. Figure [5] presents emission 
spectra calculated at different pumping rates W for large 
detuning A = ujx ~ w c = 5g, and Fig. [3] shows how the 
mean numbers of excitons (nx) and photons (nc) depend 



on pumping rate for different values of A. The figures 
reveal nontrivial behaviour of spectra and particle statis- 
tics when W changes. To understand it wc analyze the 
cigenstates of the Hamiltonian (|T|) . 

Two m-particle states can be considered for sufficiently 
large detuning (g\/m <ti A) as exciton-like (X) and cavity 
photon- like state (C). Their energies in the second order 
of perturbation theory in g are given by 



Effl = rahjje - mg 2 /A, 



(16a) 



£ T ( „ X) = (m-l)hujc+rujx+rng 2 /A, m = 1, 2 . . . . (16b) 
Only the exciton-like states 



|mX) 



1,1 



A 



m,0), 



1,2, 



(17) 



are populated by pumping to the quantum dot and they 
determine the emission spectrum. The first term in (|17[) 
corresponds to occupied dot and to— 1 photons, while the 
second decribes small admixture of the state with empty 
dot and to photons. 
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FIG. 3: (Color online) Average numbers of excitons (a) and 
numbers of photons (b) as functions of the pumping rate, cal- 
culated for the values of detuning (wx — Wc)/<? = 0, 2, 4, 5, 7, 
indicated at each curve. Other parameters are the same as 
for Fig. ffj The inset schematically illustrates the three low- 
est manifolds of Jaynes-Cummings ladder. Vertical arrows 
indicate radiative and non-radiative transitions between the 
exciton-like (X) states, governing the kinetics of the system. 



Let us consider which eigenstates determine the emis- 
sion spectra for different pumping. For very small pump- 
ing rates, corresponding to Fig. [2^., only the manifold with 



m = 1 is populated. The spectrum consists of the peak 
close to the exciton energy Ef- and very weak peak at the 
energy E^ of the almost empty photon-like state. Such 
spectrum is typical for linear in pumping regime at large 
detuning £ 

Interestingly the state with one exciton, |1X). becomes 
significantly populated already at small pumpings. In- 
deed, the pumping rate to this state is just 



ir, 



x 
o-Hix) 



= W, 



while its decay rate (determined by the photon escape 
from the cavity) is small: 



Hi 



| IX)- 



,„ = rc*(fy 



(18) 



since photonic fraction is almost negligible in this state. 
At W > WE-m-to the dot is already almost completely 
occupied while the cavity is still almost empty, (nc) ~ 
(ff/A) 2 < (n x ) ~ 1, cf. Figs.H a and b. 

By contrast, the pumping to the second state |2X) is 
inefficient, 



Tr, 



1X)->2X) 



W x 



i> « w 



while its radiative decay rate is high, 



^|2X> 



| IX) 



r c . 



Consequently, to populate second state Eq. (|17j) high 
pumping is required, contrary to the first state, see the 
inset to Fig. [3] Interestingly, that relatively high value of 
radiative rate W^xwiix) means that even for negligible 
population of the state |2X) its contribution to the emis- 
sion spectrum can be already important. The frequency 
of the corresponding peak is close to the photon mode, 
Ef - Ef = hjj c + hg 2 /A, sec Fig. [5b. This peak domi- 
nates in the spectrum when the pumping rate exceeds Tq- 
Increase of pumping rate up to 



W~T C x 



fJ 



(19) 



makes population of the second state |2X) considerable 
>2X> ~ W \2X)^\ix))- At such Pumping 



(in this case W£ x \ 



average exciton and photon numbers are both on order of 
unity, as is demonstrated by Fig.[3J The further increase of 
the pumping rate leads to the population of higher states 
(|17p with m > 3. For higher m the exciton fraction in 
states (JTTJ) is smaller so that the number of excitons slowly 
decreases and tends to 1/2, as in case of zero detuning. 
Thus, the exciton population has a plateau at (g/A) 2 < 
VF/rc < (A/g) 2 . This plateau becomes more prominent 
for higher values of detuning, see Fig. [3h. 

The ratio between the emission and pumping rates for 
m > 3 manifolds remains the same as for transition 
|2X) — > | IX), so the photon number increases linearly 
with pumping, see Fig. [5}d. The peak frequency remains 
approximately equal to ujq + 9 2 / A because the states E* 
are almost equidistant. At even higher pumping rates, 
the detuning becomes smaller as compared to the effec- 
tive exciton-photon coupling g-^/nc so that our analysis 



presented in Sec. IIII Al for high powers at zero detun- 
ing becomes applicable. Increase of pumping leads to the 
transition into lasing regime with narrowing spectrum, cf. 
Figs.Hfc and Hi. 



IV. MESOSCOPIC EFFECTS IN EMISSION 

In this Section we consider the microcavity where N > \ 
quantum dots are close to the photonic mode resonance. 
First we ignore the detuning between the excitonic and 
photonic energies and then take detuning into account. 

In case when all N quantum dots are the same (wx,i = 
wcj Tx.i = Tx, Wi = W), emission spectra at low pump- 
ing (NW <C Tq) are determined by the symmetrical, 
supcrradiant mode, where all the excitons oscillate in 
phase.— Supcrradiant mode belongs to the manifold with 
m = 1 particle, which consists of N + 1 states in total, 



-Ei,sr,± = fiwc ± y/NHg, 

-El.dark = f&C, 



(20) 
(21) 

TV — 1 of which, given by Eq. (|2ip. are dark. Thus, 
the emission spectrum has a two-peak structure with the 
splitting between the peaks enhanced by the factor \/~N. 
Increase of the pumping rate leads to the population of 
the higher manifolds with m > 2. Extra peaks arise in 
the spectrum due to the radiative transitions between the 
manifolds m = 2 and m = 1 since dark states of the 
first manifold may serve as final states for the transitions 
from the second one. If the number of dots is larger than 
unity, the transitions between second and first manifold 
become manifested at substantially lower powers as com- 
pared with a single dot case. This effect is related to the 
dark states Eq. ([2~Tj) . While the pumping rates to the su- 
perradiant and dark states are the same, the total decay 
rates -Di^dark of the dark states in the first manifold are 
much smaller since photon emission is not possible. For 
instance, if Tq is negligible, dark states decay only via 
pumping to the second manifold, W$ dark>->2 u- Conse- 
quently, they are already populated at small pumping and 
act as a reservoir for the pumping of the second manifold. 
It is similar to the role of the long- living state | IX) with 
long radiative lifetime in pumping mechanism considered 
in the previous Section, when the single dot detuned from 
the photonic mode was analyzed. 

It turns out that, the supcrradiant regime is destroyed 
by the relatively weak pumping. The emission spectra for 
N = 3 and N = 4 quantum dots presented in Figs. 0^, SJd 
at small pumping rate NW = 0.2Fc already show com- 
plex multipeak structure arising from the population of 
the higher manifolds. 

Let us now focus on the regime of large pumping. For 
W ^> Tq the typical number of cavity photons is large, 
so the electromagnetic field can be considered classicaly. 
Then the Hamiltonian ([1]) can be considerably simplified 
by replacement of the photon creation and annihilation 
operators c' and c by their expectation values. For m-th 
manifold expectation values are approximately equal to 
fm. Hence, the Hamiltonian (TTJ has the form 



H m m mftwc + \fmhgy(b i + b\) , 



(22) 
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FIG. 4: Photoluminescence spectra for N = 3 (panels a, c, e) 
and N — 4 (panels b, d, f) quantum dots for different pumping 
rates W . Panels (a) and (b), (c) and (d), (e) and (f) were calcu- 
lated for NW/Fc = 0.2, 4, 30, respectively. Other parameters 
are the same as for Fig. [Tj Energies of the transitions corre- 
sponding to the superradiant mode are indicated at panels (a) 
and (b) by vertical lines. Photon mode energy is used as the 
reference point, the spectra are nornralized to their maximum 
values. 



including only excitonic operators. To diagonalize it we 
notice, that the two-level dot behaves like the spin 1/2. 
In particular, b t + b\ can be presented as 2 si , where si 
is the operator of z projection of spin 1/2 with eigenval- 
ues ±1/2. After summation over all N dots we obtain 
the operator of total spin projection. Its eigenvalues are 
cquidistantly distributed between —N/2 and N/2, so the 
spectrum of Hamiltonian (|22|) reads 



E n 



hmujQ, + 2\i\fmhg, 



(23) 



N N N N 

l l = -7T> ~TT + !; • ■ • > IT" ~ h TT 



(a) odd N (b) even N 



Each state /i corresponds to /i + N/2 dots with "spin up" 
and —/i + N/2 dots with "spin down". State degeneracy 

equals to C^- and the total number of states is 2 N . 

The radiative transitions between manifolds m and m — 1 
conserve /i in used approximation. Decay and pumping 
rates are ^-independent. Thus, emission spectrum is de- 
termined by the states with the smallest |/i|, which have 
the highest degeneracy. It is important, that the smallest 
value of \/j,\ equals to either 1/2 or depending on the 
parity of the number of dots, so 

£m,±l/2 - £m-l,±l/2 « ^C ± 1 rj= (odd N), (24) 

Zx/Tfl 

E mfi - £ m -i,o « fkvc (even N) . (25) 

These two cases are schematically illustrated in Figs. [SJi 
and [SJd. We conclude, that for odd number of dots and 
high pumping rates [re <C NW <C W* where W* is given 
by Eq. (fT4)) ] the emission spectrum has a two-peak shape, 
like for N = 1, while of even N it consists of single peak 
at u = wc . 

Our above analysis is fully confirmed by the results of 
numerical calculation presented on Fig. |4j Panels (a)- 
(f ) show the dependence of spectra for TV = 3 and N = 
4 dots on pumping rate. As we mentioned above, the 
superradiant regime is already destroyed at relatively low 
pumping NW = 0.2Fc since the dark states increase the 
pumping efficiency to the second manifold, see panels (a), 
(b). At relatively high pumping NW — 4Tq, shown in 
panels (c) and (d), the distance between peaks for both 
N = 3 and N = 4 becomes smaller due to saturation of 
oscillator strength, while individual peaks are wider, just 
like for a single dot. However, the spectral shape strongly 
depends on the parity of N, the doublet corresponds to 
N = 3 and singlet corresponds to N = 4. Finally, at 
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FIG. 5: Scheme of the radiative transitions for odd (a) and 
even (b) numbers of quantum dots N. 



FIG. 6: (Color online) Photoluminescence spectra calculated 
taking into account inhomogeneous broadening. Thick (black) 
and thin (red) curves correspond to N — 3 and N — 4 dots, 
respectively. Solid and dashed curves were calculated for dif- 
ferent particular realizations of random spread of excitonic en- 
ergies. The value of pumping rate was NW = 4Fc- Arithmetic 
mean values of the detuning are {\u}x,i—wc\)/g ~ 0.4, 0.5, 0.6, 1 
for black solid, black dashed, red solid, red dotted curves, re- 
spectively. Other parameters are the same as for Fig. [TJ Pho- 
ton mode energy is used as the reference point, the spectra are 
normalized to their maximum values. 



high pumping rates the system is in the lasing regime, see 
Figs. 0J: and|3f. Its spectrum consists of single peak which 
width decreases with the pumping power. 

We note, that the Hamiltonian (Q]) has been already 
extensively studied for cavity with TV ^$> 1 atoms. 
Such a system is well described by the Tavis-Cummings 
modeh 22 i 23 The specific feature of the semiconductor cav- 
ity is that the number of quantum dots strongly coupled 
with photonic mode cannot be very large, TV ~ 1 . . . 10. 
Parity-sensitive mesoscopic effects demonstrated above 
become then important. 

Moreover, the inevitable disorder leading to the spread 
of the excitonic resonant frequencies cox,i makes semi- 
conductor systems particularly different from the atomic 
cavities. In order to study the disorder effect we cal- 
culated the emission spectra of three and four quantum 
dots in a microcavity for the different realizations of the 
exciton resonance frequencies, see Fig. [6] For all four 
curves in Fig. [5] the mean arithmetic values of the de- 
tuning (l/TV)^ i=1 \uJx,i — wc| were in the order of the 
light-matter coupling constant g. Although the spectra 
are modified by disorder, the doublet (TV = 4) can be still 
well distinguished from the singlet (TV = 3). Thus, we 
conclude that the parity-dependent mesoscopic effects in 
emission spectra are relatively stable against the disorder. 



V. CONCLUSIONS 

To summarize, we have developed a kinetic theory of 
the non-linear emission of quantum dots embedded into 
the semiconductor microcavity. Wc considered the case 
of the strong coupling where the light-matter interaction 
constant is larger than the damping rates of the cavity 
and exciton, governed by the photon escape through the 
mirrors and exciton non-radiative decay, respectively. In 
this case the eigenstates of the system belonging to the 
manifolds with the total particle number m are well de- 
fined so that the emission spectra can be found from 
the Fermi golden rule. The populations of the coupled 
photon-exciton states are determined by the kinetic equa- 
tion which takes into account both pumping and decay 
processes. The linewidths of the emission spectra are 
found as corrections to the kinetic equations. 

Our method was tested in the case of a single quantum 
dot in a microcavity being resonant with the photonic 
mode. In the strong coupling regime we reproduced the 
results of Ref. \Q The transition from the linear regime 
with Rabi doublet via multipeak spectrum to the nar- 
row emission line being the characteristic of the laser is 
demonstrated. 

We have studied in detail the case of the single dot de- 
tuned from the cavity mode. At low pumping rates the 
emission is dominated by the quantum dot line. An in- 
crease of the pumping results in the population of the sec- 
ond manifold and the emission line is close to the cavity 
position. Even higher pumping results in the line narrow- 
ing and transition to the lasing regime similarly to the 
case of the resonant quantum dot. 

The cigenmodes of the interacting cavity-quantum dots 
system can be found analytically in the case of large pho- 



ton number where the cavity field is classical and the light- 
matter coupling can be formulated in the terms of angu- 
lar momentum. In the situation where several dots are 
embedded into the microcavity the emission spectra are 
shown to demonstrate parity effect: at relatively strong 
pumping there are one or two peaks in the spectrum de- 
pending on whether the dot number is even or odd, re- 
spectively. This parity effect is relatively robust to the 
spread of exciton resonance energies caused, e.g., by the 
disorder. 

To conclude, the non-linear emission spectra of the 
quantum dots in a cavity provide important information 
both about light-matter interaction and about the quan- 
tum dot ensemble. 
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Appendix A: Analytical solution for single QD in 
microcavity 

In this Appendix we present analytical solutions for po- 
lariton distribution function and emission spectrum of sin- 
gle quantum dot, strongly coupled with microcavity pho- 
ton mode. The detuning is neglected, loq = lux- 

We consider relatively high pumping, W 3> Tq , so that 
the spectrum is determined by high manifolds, m S> 1. 
Moreover, the strong coupling regime is assumed, i.e. g ~3> 
Tq, y/mg ~> W, the exciton decay is neglected. In this 
case only states with the same symmetry are coupled by 
the photon annihilation operator: 



^m—l,ii;mfi' 



(Al) 



Clearly, under these assumptions states /i = + and 
fj, = — are equally populated. First we need to deter- 
mine the stationary distribution function of polaritons 
fm = fm,+ + fm- = 2/ r „ : + = 2/ m _. Under the as- 
sumption (|A1[) kinetic Eq. (0} reads 

— j- = -T c mf m + r c (m + l)/ m +i + -yC/m-i - fm)- 

. (A2) 
Since typical values of m are high, the discrete function f m 

can be replaced by continuous distribution f(m). Finite- 
difference equation (|A2[) then, at steady state, becomes 
differential equation 

( m +( m ))^4 + 2(m-(m))^ + 2/ = 0, (A3) 



dm 2 



dm 



where (m) = W/{2T C ). The solution of Eq. (|A3]) . describ- 
ing the distribution function, reads 



/(m)oc (m+(m)) 4 < m > +1 e- 2m 



(A4) 



where the normalization constant to be determined from 
Eq. ([5]) is omitted. Thus, for high pumping, when 
(to) ^> 1, the distrubution function has Gaussian shape, 



f(m) cx exp(^- 



(m — (m))' 

2M 



(A5) 



with (to) being the average number of polaritons. 

Now we proceed to the calculation of emission spec- 
trum. In the strong coupling regime the only relevant 
off-diagonal components of the density matrix are those 
with /i = /i', which can be presented as g m ,± ;m -i,± = 
(a m ± ib m ) cxp(— iujet) , where a m and b m are real coeffi- 
cients. Eq. §5§ leads to 

W 
a m (t) = g m b m (t) + — [a m -i(t) - a m (t)\ + {T c a) m , 

W 
b m (t) = -g m a m {t) - —b m (t) + {T c b) m , (A6) 



where g m = g(y/rn — \/m — 1) and the term 



(r c d) m = r c 



y/m(m + Tjdm+i - (m - -)d Tl 



,d = a,b, 

(A7) 



describes the cavity decay. Initial conditions (fTT|) are 

a m (0) = Vmf m , b m (0) = (A8) 

and emission spectrum (fTO)) is given by 

(A9) 



I(u>) cx 2, \fin / cos[(w — u!c)t]a m (t)dt 

_ Jo 



Similarly to kinetic equation, system (|A6|) can be solved 
efficiently by introducing continuous functions a{m), 6(m), 
which satisfy differential equations 



da Px da da 



(A10) 



dt y 2 



Keeping in mind initial conditions 
tions of (jAlOp in the form 



dm 

_dfc 

dm 
we seek the solu- 



a(m) = y/mf(m)A(t), b(m) = y/rnf{m)B(t) . (All) 



To proceed further we integrate (|A10[) over to using iden- 
tities 



/ 



dma(m) m 
dmm 3/2 



\/{m)A, I dmy/ma(m) sa (m)A , (A12) 
da 3(m) f , , — da A 



3(m) f — da 

— A J dmVm-— 

A 



2' 



dm(y/rn — y/m — l)a(m) ~ , 



which follow from Eqs. (|A5|) . (|A11[) in the leading order 
over (to). The functions A(i), B(t) satisfy the coupled 
differential equations 



dA _ dB 

d^ = • 9offB ' W 



-g eS A - 2jB 



(A13) 



with initial conditions A(0) — 1, B(0) = 0, where g e ff 
g/(2y / (m)) and 7 = (W + r c )/4. The solution reads 



A(t) = exp(— 7*) ( cos Oi + -^ sin Oi 



(A14) 



with O = V^S ~~ T 2 • Finally, the emission spec- 
trum is obtained by Fourier transformation, /(w) cx 
(m) J cos[(w — ujc)t]A(t)dt . 

At a threshold value of the pumping, W — W* = 
2g 2 / 3 r c the Q turns to zero. The spectrum has differ- 
ent form, depending on the relation between the pumping 
rate W and W* . For W < W* the value of SI is real and 
the spectrum has two-peak shape 



I(w) oc 



2fi- 



•wc 



2n 



WC 



(w - Wc - ^) 2 + 7 2 (w - w c + n) 2 + 7 2 ' 

(A15) 
The peak width is increasing with pumping. At relatively 
low (r c < W < W*) pumping rates (|A15[) can be sim- 
plified to Eq. (flU|) . Above the threshold the spectrum has 
single peak shape 

1 1 

/(w) oc 



(w-wc) 2 + (|0|- 7 )2 (w-w c ) 2 + (|f!|+7)2- 

(A16) 
Only the first term in Eq. (|A15[) is of interest at high 
pumping {W ^S> W*), leading to the peak (TTSl) . narrowing 
when W increases. We remind that our theory is valid at 
W < ,g 2 /r c ~ W* x (g/r c ) 4/3 . At very high pumping 
rate the laser self-quenching takes place.^21 



Appendix B: Effects of the blueshift, random sources 
approach 

In this Appendix we consider the effects of exciton- 
exciton interactions on the emission spectra of quantum 
dots embedded in the microcavities. In order to elucidate 
the role of interactions we use classical model of Ref. |fj and 
treat excitons and photons as classical oscillators (quan- 
tum treatment for the model case of interactions is pre- 
sented in Ref. [24| and for the lasing regime in Ref. |25I ). 
These oscillators are conveniently described by the di- 
mensionless exciton polarizations Pi = (bi) (i = 1, . . . ,n) 
and dimensionless electric field amplitude in the cavity 
E = (c). Corresponding linear equations of motion can be 
derived from Hamiltonian Eq. ([T]) by considering Heisen- 
berg equations of motion for the excitonic and photonic 
operators, taking average of the Heisenberg equations and 
neglecting high order correlators of exciton fields. So 
far, exciton-exciton interactions were disregarded. In the 
scmiclassical approach they can be introduced as non- 
linear terms in the coupled oscillators model: 



dE 
~dt 



iw c 



r c 



E + gJ^Pi, 



dt 



1WA". 



r x 



fl + ^+E^i^i 2 \Pi 



gE + Wi(t), i = l,...,N. 



(Bla) 



(Bib) 
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Equation for the cavity field remains linear while the equa- 
tions for the polarizations, Eq. (|Blb[) . acquire anharmonic 
contributions described by the real coefficients oiji which 
determine energy shift of ith exciton due to the interac- 
tion with jth exciton. Rigorous derivation of P 3 terms 
and discussion of their microscopic origin is out of the 
scope of this paper. Terms Wi(t) in Eq. (jBlbj) are the ran- 
dom forces which determine exciton generation in quan- 
tum dots£ 

In order to make the analysis more transparent we, in- 
stead of studying Eqs. (|B1|) in their complexity focus on 
the single exciton in one quantum dot and take into ac- 
count simplest possible interaction in the form a|P| 2 P. 
The frequency spectrum of such an oscillator driven by 
the white-noise random force 

(w(t)w{t')) =WS(t-t r ), 

where the brackets denote averaging over the different re- 
alizations of the pumping, is related to the joint probabil- 
ity T{P, t;P', £'), which can be rigorously determined by 
Fokker-Planck equation tectmique i 26 i 27 

Similarly to the Eq. (|9]) in the quantum case, the spec- 
trum reads 

/>oo 

(|PH| 2 ) oc Re / dt (P(i)P*(0)) e iw * (B2) 

Jo 

where the correlator (P(i)P*(0)) is given by 



(P(i)P*(0)) 



PP'*J r (P 7 t;P',0)d 2 Pd 2 P', (B3) 



and d 2 P = dRePdlmP. Joint probability can be ex- 
pressed as^ 6 - 



P{P, t; P', 0) = g(P, t\P', O)F(P') 



(B4) 



via the conditional probability Q(P,t\P' ,0) and 
the stationary distribution function -P(P') = 
\im t ^ooG(P',t\P",0). The conditional probability 

satisfies Fokker-Planck equation 26 



dg(p,t\p',o) 

dt 



£(P,P*)£(P,£|P',0), (B5) 



g(p,o\p',o) = s(P-p'), 



C(P,P*) 



i( L o a +a\P\ 2 ) + ^- 



d a 2 

P + cc +W 

dP ' + + dPdP* 



allows analytical solution for arbitrary 
— For relatively high pumping rate the spectrum 



Eq. 

pumping 

is nonzero only for for u> 



wx and takes the form 



(|P(w)| 2 ) (xW(oj-oj x )cxp 



LU — U>X 

a(n x ) 



, u > ujx- (B6) 



We introduce the parameter (nx) = W/Tx which char- 
acterizes the stationary excitonic population (JP(£)| 2 ), it 
depends both on pumping strength and on the excitonic 
decay rate, but it does not depend on the value of the 
non-linearity. It follows from Eq. (|B6|) that the emission 
spectrum of the system is strongly asymmetric, the maxi- 
mum position is shifted by a (nx) from the non-interacting 
position. Moreover, the spectrum is strongly broadened: 
its width is of the same order as the energy shift. Eq. (|B6|) 
can be understood as follows. For the strong non-linearity, 
where a (nx) ~S> Tx the shape of the spectrum is de- 
termined by the fluctuations of the excitonic population 
nx~ The steady state distribution of nx is the same as 
for a = 0, i.e. it is described by the Poisson formula. In 
these case the root mean square value of fluctuations of 
nx is proportional to the average population (nx)- The 
spectrum is also asymmetric since interactions blue-shift 
exciton energy only, hence, high energy wing is larger than 
low energy one. 

The very same considerations can be applied for the sys- 
tem of coupled dots and microcavity. The pumping results 
in the blue shifts of exciton energies and increase of the ex- 
citon resonance widths. As a result, the light-matter cou- 
pling becomes weaker. It is worth to stress, that the blue 
shift and broadening are contributed by all particles in 
the system, including the excitons generated in the higher 
energy states and in the wetting layer. For rather realistic 
pumping rates the blue shift can reach about 1 meV (see, 
e.g., Ref. [29 ). hence the broadening of the exciton state 
due to the particle number fluctuations can exceed the 
Rabi slitting being about 0.5 meV in the state-of-the-art 
structures^ Therefore, the observation of the non-linear 
effects for quantum dots embedded into the microcavitics 
may be hindered by the particle number fluctuations. 
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